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Abstract 

This paper proposes a form of MPC in which the control variables are moved 
asynchronously. This contrasts with most MIMO control schemes, which as- 
sume that all variables are updated simultaneously. MPC outperforms other 
control strategies through its ability to deal with constraints. This requires on- 
line optimization, hence computational complexity can become an issue when 
applying MPC to complex systems with fast response times. The multiplexed 
MPC scheme described in this paper solves the MPC problem for each sub- 
system sequentially, and updates subsystem controls as soon as the solution is 
available, thus distributing the control moves over a complete update cycle. The 
resulting computational speed-up allows faster response to disturbances, which 
may result in improved performance, despite finding sub-optimal solutions to 
the original problem. 
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1 Introduction 



1.1 The basic idea 

Model Predictive Control (MPC) has become an established control technology in 
the petrochemical industry, and its use is currently being pioneered in an increasingly 
wide range of process industries [231 El]- It is also being proposed for a range 
of higher bandwidth applications, such as ships [22], aerospace [211 [2l|, and road 
vehicles |20j . This paper is concerned with facilitating applications of MPC in which 
computational complexity, in particular computation time, is likely to be an issue. 
One can foresee that applications to embedded systems, with the MPC algorithm 
implemented in a chip or an FPGA [3 [131 [U], are likely to run up against this 
problem. 

MPC operates by solving an optimization problem on-line, in real time, to determine 
a plan for future operation. Only an initial portion of that plan is implemented, and 
the process is repeated, re-planning when new information becomes available. Since 
numerical optimization naturally handles hard constraints, MPC offers good per- 
formance while operating close to constraint boundaries [18]. Solving a numerical 
optimization can be a complex problem, and for situations in which computation is 
limited, the time to find the solution can be the limiting factor in the choice of the 
update interval. Most MPC theory to date, and as far as we know all implementa- 
tions, assumes that all the control inputs are updated at the same instant. Suppose 
that a given MPC control problem can be solved in not less than T seconds, so that 
the smallest possible update interval is T. The computational complexity of typical 
MPC problems, including time requirements, tends to vary as 0((m x Ny)'^), where 
m is the number of control inputs and is the horizon length. We propose to 
use MPC to update only one control variable at a time, but to exploit the reduced 
complexity to update successive inputs at intervals smaller than T, typically T/m. 
After m updates a fresh cycle of updates begins, so that each whole cycle of updates 
repeats with cycle time T. We call this scheme multiplexed MPC, or MMPC. We 
assume that fresh measurements of the plant state are available at these reduced 
update intervals T/m. The main motivation for this scheme is the belief that in 
many cases the approximation involved in updating only one input at a time will 
be outweighed — as regards performance benefits — by the more rapid response 
to disturbances, which this scheme makes possible. It is often the case that "do 
something sooner" leads to better control than "do the optimal thing later". Fig. 
[T] shows the pattern of input moves in the MMPC scheme with m = 3, compared 
with the conventional scheme in which the three input moves are synchronized. We 
will refer to conventional MPC as Synchronized MPC, or SMPC, in the rest of this 
paper. 

The scheme which we investigate here is close to common industrial practice in 
complex plants, where it is often impossible to update all the control inputs simul- 
taneously, because of their sheer number, and the limitations of the communications 
channels between the controller and the actuators. 

In addition to treating the 'nominal' MMPC case, in which the model is assumed to 
represent the plant perfectly, we extend MMPC to guarantee robust constraint satis- 
faction and feasibility of all optimizations despite the action of unknown but bounded 
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Figure 1: Patterns of input moves for conventional 'synchronized' MFC (left), and 
for the Multiplexed MFC (right) introduced in this paper. 



disturbances. These are key issues in MFC: performance benefits are achieved by 
operating close to constraint boundaries, but when the state evolution no longer 
matches the predictions, constraint violation and infeasibility can result. Many 
methods have been developed to endow conventional synchronous MFC with ro- 
bustness [m [1]. For use with MMFC, we have adopted the constraint tightening 
approach |1H [HI [271 [25] , in which the constraints of the optimization are modified to 
retain a margin for future feedback action. Since only the constraint limits are mod- 
ified, the computational complexity remains the same as for the equivalent nominal 
MFC. Constraint tightening is therefore well-suited to MMFC, which is aimed at 
computation- limited applications. 

Various generalizations of our scheme are possible. For example, subsets of control 
inputs might be updated simultaneously, perhaps all the inputs in each subset being 
associated with one subunit. The assumption of equal intervals between the updates 
of plant inputs is not essential to the MMFC idea. Any pattern of update intervals 
can be supported, providing that it repeats in subsequent update cycles. A further 
generalization, albeit involving a significantly harder problem, would be not to up- 
date each control input in a fixed sequence, but to decide in real time which input 
(if any) needs updating most urgently — one could call this just-in-time MFC. 

1.2 Related Work 

MMFC is related to distributed MFC (DMFC) [7], both dividing the optimisation 
into smaller sub-problems. Several works have been published which propose 'dis- 
tributed MFC in the sense that subsets of control inputs are updated by means of 
an MFC algorithm. But these usually assume that several sets of such computa- 
tions are performed in parallel, on the basis of local measurements only, and that 
all the control inputs are then updated simultaneously. In some applications, such 
as formation flying of unmanned vehicles [lO], it is assumed that the state vectors 
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of subunits (vehicles) are distinct, and that coupHng between subunits occurs only 
through constraints and performance measures. In [33] five different MPC-based 
schemes are proposed, of which four are distributed or decentralized MPC schemes 
of some kind. Their schemes 4 and 5 are the closest to our multiplexed scheme. In 
these schemes an MPC solution is solved iteratively for each control input, but it is 
assumed that no new sensor information arrives during the iteration, and that all the 
control inputs are updated simultaneously when the iterations have been completed. 

Various approaches to robust DMPC have been investigated, including worst-case 
predictions [12], retention of "emergency" plans [311 |6], invariant "tube" predic- 
tions |32| and constraint tightening |26] . Work on DMPC has typically focussed 
on spatially distributed systems with some structure in the system, e.g. teams of 
vehicles with decoupled dynamics. In contrast, our new robust MMPC makes no 
assumptions on the overall system structure, and considers temporal distribution, 
breaking the optimisation down into a sequence of smaller problems, potentially on 
the same processor. 

In [2] a similar scheme to ours is proposed, but it is assumed that a limitation occurs 
on network bandwidth rather than on central computing resources. Hence optimal 
trajectories are computed for all the plant inputs, but these are communicated to the 
plant one input (or one group of inputs) at a time — with the optimization taking this 
communication restriction into account. If the communication sequence is fixed and 
periodic then the scheme proposed in |2j is essentially the same as a version of MMPC 
to which we previously referred as 'scheme 1', except that we allowed constraints on 
inputs and states [16]. |2] also considers the case that a feedback law is fixed for each 
input (or group of inputs), that the inputs are updated according to some periodic 
scheme, and that a heuristic is used to determine (online) the best point in the period 
for a given state; this gives a heuristic version of 'just-in-time control' as defined 
above, though not really MPC any longer, since the feedback law is assumed to 
be predetermined. We emphasize that the driving factor behind the development of 
MMPC is operation in a processor- limited environment, motivating decomposition of 
the optimisation to reduce computational delays. Therefore we have not considered 
the impact of communication limits: indeed, in many applications of MMPC, the 
computation may take place serially on a single processor, and thus communication 
is not a concern. 

In |3D] MPC is considered with opposite assumptions to ours on update rates. There 
the plant inputs are considered to be updated relatively frequently, compared with 
the rate at which output measurements become available. This is in contrast to 
MMPC, in which the plant outputs are assumed to be measured relatively fre- 
quently, compared with the rate at which inputs are updated. It is remarked in |30] 
that the predictive control law which results (with the specific assumptions made 
there) is periodic, the period being the ratio of the input update rate to the output 
measurement rate (assuming this is an integer). A similar observation is central to 
the development in our section [3j 

An alternative strategy for speeding up the computations involved in MPC is 'ex- 
plicit MPC, which involves off-line precomputation of the 'pieces' of the piecewise- 
affine controller which is the optimal solution |20j. But that is not feasible if the 
number of 'pieces' required is excessively large, or if the constraints or the plant 
model change relatively frequently. 
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MMPC was introduced by us in [16]. Robust MMPC was first described in ^8]. In 
p9j our MMPC idea was applied (by others) to the control of an aircraft engine. In 
|15| an experimental evaluation of MMPC is reported. 

1.3 Structure of the paper 

The rest of this paper is organized as follows. In Section [2] a formulation of MMPC 
is presented in detail. Section [3] establishes the nominal stability of MMPC with this 
formulation. Section |4] then derives a formula for the value of the cost function at- 
tained by MMPC. Section [5] develops a more elaborate formulation of MMPC, with 
the objective of guaranteeing robust feasibility, and establishes an appropriate the- 
orem. Section [6] gives numerical simulation examples and compare the performance 
of MMPC with SMPC for cases with significant plant uncertainty, represented by 
unknown but bounded disturbances. Finally, concluding remarks are given in Sec- 
tion [7l 



2 Problem formulation 
2.1 Preliminary 

We consider the following discrete-time linear plant model in state-space form, with 
state vector Xk € M" and m (scalar) inputs ui^fc, . . . , Um,k' 

m 

Xk+1 = Axk + ^ BjAuj^k (1) 
j=i 

where each Bj is a column vector and Auj^k = uj^k — Uj,k-i- (This could be general- 
ized to the case where Bj € M^^W and Auj^k £ , with Y2jPj inputs.) We assume 
that (A, [Bi, . . . , -Bm]) is stabilizable. For ease of notation, when we drop the index 
j, we mean the complete B matrix and the input vector so that the system ([T]) may 
be written as 

Xk+i = Axk + BAuk 

We assume that at time step k the complete state vector Xk is known exactly from 
measurements. We will consider only the regulation problem in detail, but tracking 
problems, especially those with non-zero constant references, can be easily trans- 
formed into equivalent regulation problems [3, sec. 3. 3]. 

Multiplexed MPC, at discrete-time index k, changes only plant input Au„^k),ki where 
cr(k) is an indexing function which identifies the input channel to be moved at each 
step, and is defined as: 

a{k) = {k mod m) + 1 (2) 

(We assume, without loss of generality, that we update input 1 at time index 0.) 
The asynchronous nature of the multiplexed control moves, as illustrated in Fig. [T| 
is captured by the constraint 

Auj^k = if j / a{k). (3) 
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It is then possible to rewrite the system dynamics Q as a hnear periodically time- 
varying single-input system: 



Xk+i = Axk + -B^(A,.)Attfc 



(4) 



where Auk = Au, 



a{k),k- 



From this point onwards, we use this periodic description of 



the plant so that we can draw on known results for periodic time-varying systems. 



Remark 1 Some of the generalizations to which we alluded in section \Ll\ could he 
treated by redefining the sequencing function a(-) appropriately. For example, for a 
particular 3-input system, updating the inputs in the sequence (1, 2, 1, 3, 1, 2, 1,3,.. .), 
thus updating one of the inputs twice as often as the others, could be represented in 
this way. 



The unique advantage of MPC, compared with other control strategies, is its capacity 
to take account of constraints in a systematic manner. As usual in MPC, we will 



suppose that constraints may exist on the input moves, Auk € 



Jcr(fc) 



and on states. 



Xk € X, where X and ^^(^k) ^-^^ compact polyhedral sets containing the origin in their 
interior. Note that the control move set depends on the time, since the channel to 
be moved differs from step to step. If constraints on the actual control inputs u are 
required, then u must appear in the augmented state x, and those constraints can 
be incorporated in the state constraint set X. 

Let = {Nu — l)m + 1 where is the control horizon, a design parameter which 
will later be used to denote the number of control moves to be optimized per input 
channel of the original system ([T]). The A^-step prediction model at time k for the 
system described by Q is 



Xk+l\k = ^Xk\k + Ga{k)^Uk\k 



(5) 



where 



^k+l\k 



Xk+l\k 
Xk+2\k 

Xk+N\k 



AU, 



k\k 



^Uk\k 
A'Ufc+llfc 

^Uk+N-l\k 





A 















Ga{k) — 



B 



AB. 



•a{k) 



a{k) 



B. 





'<T(k+l) 








A^'-^B 



(6) 



^^'fc+ilfc denotes the prediction made at time /c of a control move to be executed at 
time k + i, and denotes the corresponding prediction of x^+j, made at time k. 



2.2 The MMPC Algorithm 

In the following, -fCo-(fc) denotes a pre-specified stabilizing linear periodic state- 
feedback controller of Q; {^Xi{K„(^i^-^)) denotes a sequence of sets in which none 
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of the constraints is active, and which satisfies the 'periodic invariance' condition 
for the linear periodic system Q when the feedback controller 

Auk = -K„(k)Xk (7) 
is applied, namely 

Xk G Xi{K„(^k)) -K„^i,)Xk G U^(fc) and [A - B^(^i,)K„(^k))xk G 
and of course Xj{Kf,.(^f^)) C X, for cr{k) = 1, . . . ,m. 

Some assumptions must be made about those inputs which have already been 
planned but which have not yet been executed. We will assume that all such planned 
decisions are known to the controller, and that it assumes that they will be executed 
as planned, i.e., 

^Uk+i\k = Auk^i\k-i, z/0,m,2m, ... (8) 

(In fact, new decisions will be made at time k + i in the light of new measurements.) 
MMPC solves the following finite-time constrained linear periodic control problem: 



Jk — Fcr{k){Xk+N\k) + Silo^ (lkfc+j|fcllg + II '^'"fc+j|A; llr ) 

^Uk+i\k, {i = 0,m,2m,...,N -1) 
Aufc+i|fc G U^(fc+i), {i = 0,...,N -1) 
Xfc+j|fcGX, (i = l,...,iV- 1) 
Xk+N\k G Xi{K^(^k)) 
Xk+i+i\k = + i?o-(fc+i)An,fc+j|fc 

^Uk+i\k = Atifc+i|fc_i, (i / 0, m, . . . , iV - 1) 

(9) 

where -F(T(fc)(^fc+Ar|fc) > is a suitably chosen terminal cost. 

We denote the resulting optimizing control sequence as Au°(xk)- Only the first 
control Au^ in Au°{xk) is applied to the system at time k, so that we apply the 
predictive control in the usual receding-horizon manner. 

In MMPC, there are essentially m MFC controllers, operating in sequence, in a 
cyclic manner. They share information, however, in the sense that the complete 
plant state is available to each controller — although not at the same times — and 
the currently planned future moves of each controller are also available to all the 
others. 

For clarity, we set out the following algorithm which defines 'nominal' MMPC (as 
contrasted with 'robust' MMFC which will be introduced in section [s]): 

Algorithm 1 (Nominal MMPC) 

1. Set k := kf). Initialise by solving problem but optimising over all the 
variables An;j_|_j|fc, i = 0, 1, . . . , — 1. 

2. Apply control move Au^(^f.-j^i^ = Au^j^ 

3. Store planned moves A-u^ ^^i^. 



'P<T{k){xk) : Minimise 
wrt 
s.t. 
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4- Pause for one time step, increment k, obtain new measurement x^. 

5. Solve problem Q. 

6. Go to step\^ 

Note that Step [T] involves solving for inputs across all channels, not just channel a{k). 
This type of initialisation requirement is common in distributed MPC. Subsequent 
results do not depend on the optimality of this initial solution, only its feasibility. 



3 Stability of MMPC 

In this section we establish sufficient conditions under which the MMPC scheme gives 
closed-loop stability. We then apply standard results on optimal control of periodic 
systems to our plant written in the form of Q, assuming that all constraints are 
inactive, to propose a terminal cost F{-) which, when used in the MMPC algorithm 



introduced in section 2.2 ensures stability of the closed loop even when constraints 
are active. 

Theorem 1 MMPC, obtained by implementing the nominal MMPC Algorithm [7], 
gives closed-loop stability if the problems are well-posed, and if the set of terminal 
costs {Ff^{-)} satisfies 

F^+{[A- B^+K^+]x) + \\x\\l + \\K^+x\\f. < F^{x) for a = l,...,m. (10) 

where = {a mod m) + 1, namely the cyclical successor value to a. 

Proof: 

The proof follows a standard argument for MPC stability proofs (see for ex- 
ample), adapted to our setting. It depends on the constrained optimization being 
feasible at each step, and the feasibility at any particular time step depends on the 
details of the constrained optimization problem that is being solved. For the nomi- 
nal case with a perfect model and in the absence of disturbances, if feasible solutions 
are obtained over an initial period, then feasibility is assured thereafter. 

Let An^_|_^|^ denote the optimal solution to ^ at time step k, for i = 0,m, 2m, . . . , A^— 
1, let denote the corresponding state sequence, for i = 1, 2, . . . , A^, and let 

be the corresponding value of the cost function J^,. 

Then a candidate input sequence to be applied to the plant at the next time step, 
-|- 1, is 



^%+l+m|A;+l-m' '^%+2+m|A;+2-m' ' ' • ' ^%+2m|fc' • • ■ ' 
^'^k+l+{Nu-2)m\k+l-m^ ^%+2+(Af„-2)m|A;+2-m' • • • ' 

^^k+{N^-l)m\k^ ~^ff(fc+l)^fc+Af|fc) (11) 

(recall that — 1 = (Nu — l)m). The input sequence applied at time step k is the 
same, except that the initial term Ail^|^ is pre-pended to it, and that the final term 
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I^a{k+i)^*k^N\k omitted. Let the cost obtained with the candidate solution (11) 
be Jk+i- Then 



"^fe+l - Jk - ll^fT(fe+l)2;fc+Ar|fellr + -^CT(fc+l) \ {^- Bcr{k+l)Ka(k+l))xlj^N\k^ 

~ ll^^fc|A:llr - ~ FcT(k){Xk+N\k) + \\Xk+N\k\\1 (12) 



Hence Jk+i — Jk < if (10) holds. Now optimisation at time step k + 1 will result 
in a value function 

Jk+i < Jk+i (13) 
and hence 

JUi < J*k (14) 



if (10) holds. 



But Jfc > for aU k, hence J^^^ - ^ 0. But, from ([12|)-([14|) we have that 



Jk+l Jk — \\^*k\k\\q ll^'"fc|fellr (15) 

Hence (and Au*|^ 0). But Xk\k = Xk, so Xk 0. 



Remark 2 Note the implicit assumption that N is chosen sufficiently large to en- 
sure feasibility of the constrained optimisation problem posed. Also note the assump- 
tion in each planning optimisation that the linear state feedback law ([T]) is applied 
at every step after the end of the optimisation horizon. 

The following results on unconstrained infinite-time linear quadratic control of peri- 
odic systems are known j4j. Consider the plant (|4| and the quadratic cost function 



«^fe = X] (lkfc+i||g + ||Anfc+,||2) (16) 

i=0 

Then this cost is minimised by finding Pi, i = 1, . . . ,m, the Symmetric, Periodic and 
Positive Semidefinite (SPPS) solution of the following discrete-time periodic Riccati 
equation (DPRE) 

Pk = A^Pk+iA - A^Pk+iB^^k){Bj^k)Pk+iB^ik) + r)-isJ(,)Pfc+iA + q (17) 
and setting 

Auk = -K„(k)Xk (18) 

where 

Ka{k) = (^J(fc)-Pf7(fc+i)^^(fc) + r)"^5j(^)P^(fc+i)^ (19) 
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Furthermore, the resulting minimal value of Jk is given by = x|'Po-(A:)2;fc- Thus 
one way of choosing a suitable set of terminal costs to satisfy (10) is to set 



which leads to 



J, 



fe+i 



T* 



\X 



k\k\\q 



lAti 



(20) 



(21) 



Remark 3 The terminal cost (20) would he the optimal cost-to-go if at each step 



k the optimisation was over future values of all input channels, rather than those 
in input channel (j{k) only. A version of MMPC in which this is done was called 
'scheme 1 MMPC in our earlier paper Ji6]/ . The version presented in this paper was 
called 'scheme 2' in flW . We no longer advocate 'scheme 1\ as it does not give any 
reduction of computational complexity, compared with conventional SMPC. 



4 Cost of MMPC when constraints are inactive 



Each solution of the optimisation problem ^ depends on the plans made in previous 
optimisations. Hence the optimal cost obtained with MMPC, even in the case that 



all constraints are inactive, is not given by (20). In this section we will introduce 



an augmented state which includes those existing plans that are not going to be 
modified by the current optimisation. This will allow us to obtain, in Theorem 
[2| an expression for the optimal cost of the same form as (20). This will provide 
an analysis tool for predicting and comparing the performance of various MMPC 
designs. In the process we will see that MMPC can be rewritten in a more familiar 
MPC form, but with a periodically time-varying (augmented state) model. 

Note that a similar development could be used to compute the optimal MMPC 
cost if the set of active constraints was constant and known. The nature of the 
MMPC control law in such circumstances is also linear periodic. Consequently the 
MMPC control law in the presence of constraints is piecewise-linear-periodic; as in 
the standard 'explicit' MPC case, the 'pieces' correspond to regions of the state 
space in which the set of active constraints remains constant. 

The development of this section, in particular Theorem [2j facilitates performance 
evaluation of MMPC in certain circumstances. For example, it is useful for evaluat- 
ing the trade-off between the restricted optimisation performed by MMPC and the 
reduced update rate available with conventional SMPC. 



4.1 Unconstrained MMPC as periodic state feedback 

We introduce the following definitions, which gather together those variables which 
are optimised at each step by the MMPC algorithm: 



'^Uk,i\k 



^^k+i\k 

^"^k+m+ilk 
^Uk+{Nu-2)m+i\k 



(22) 
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for i = 1,2, ... ,m. 



Au 



k+(Nu-l)m\k 



(23) 



Recall that N = (N^ — l)m + 1 where A^^^ is the control horizon, a design parameter 
which denotes the number of control moves to be optimized per input channel of the 
original system (IT). By grouping the predicted control signals into m vectors, the 
prediction model (51) can be re- written as 



\k 



(24) 



where Auf^^^j^ and Ano|fc are as defined in (22) and (23), respectively, and gf^'^^ (i = 
1, . . . , m) are matrices whose columns are columns of the Go-(fc) matrix ([6]), namely. 



9i 



'^^^^ is the matrix whose columns are columns 1, 1 + m, . . . , 1 + {Nu — \)m of the 



matrix Go-(fc) > while g^^^^ ,{i = 2, . . . ,m) contains columns i,i + m, . . . ,i + {Nu — 2)m 



columns of Gcr(k)- 

In MMPC only Auq^^ is taken as the decision variable at time k, and appropriate 
assumptions are made about Auk^i\k, i = 1, . . . , m — 1. Note that the length of AuQ^k 
is Nu while the length of Auk^i\k for i = 1, . . . , m — 1 is Nu — 1. When Nu = 1, 
^''^k,i\kj i = 1, ■ . ■ ,Tn — 1, become zero-length vectors. In MMPC we assume that 
the Auk^i\k, i = 1, . . . , m — 1 are those inputs which have already been planned in 
previous steps but have not yet been executed, namely 



'^Uk,i\k = ^^A:-l,i+l|fc-l> {i 



1, 



, m 



!)• 



(25) 



We define the vector Au^^j^_^ which holds the previously planned but not yet exe- 
cuted control moves as 



^^fc-l,3|fc-l 



Auk. 



l,m\k-l 



(26) 



Thus, it can be deduced from (|24|) that, if no constraints are active, then the MMPC 

(27) 



control law is a linear periodic state feedback: 
where we have introduced the augmented state vector 
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Xk 

^<lk-. 



(28) 
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4.2 A formula for the MMPC Cost 



Using the augmented state vector introduced in (28), the dynamics of the plant 
operating under MMPC can be expressed as 



(29) 



where 



A 



A 
Au 

1 
0/ 



B. 





Ba{k) 






B^ 




" 




' " 












Bu = 









/ 



The value of the quadratic cost (see also Remark |4]) 

oo 



(30) 



i=0 



when control law of the form (27) is applied, is given by the following theorem: 



Theorem 2 The value of the cost (30) for the system (29), when any stabilising 
linear periodic state feedback of the form (27) is applied, is given by 

Jk = Ji,k + Ju,k = S,k P^+u,a{k)S,k (31) 

where J^,^ = ilP^^a{k)ik and Ju,k = ilPu,a{k)ik, and Ps^^aik), Pu,a{k) and P^+u,a(k) 
are, respectively, solutions of the following Lyapunov equations: 



\.a{k) - *a(fc)^C,'T(/c)*(T(/c) + *J(fe)'9*<T(fc) 

T{k) = *a(fc)A«,0-(fe)*£T(fc) +^J(fe)^^(T(fe) 



+^l(k)PK„(k) 



(32) 
(33) 

(34) 



where 



a{k) 
aik) 



■ ^a{k+l)^a{k) 



a(fc) 



aik) 

Q 

R 



a{k} 
^a(k+l)^a(k} 



^cy{k+m~l) ■ ■ ■ ^cr{k) 

k, 



aik) 
Ka(k+l)^a(k) 



a{k+m-l)^a{k+m-2) 



a{k) J 



diag{q, q,--- ,q) 
diag{r, f, • • • ,f) 



ana 
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Proof: 

The closed- loop dynamics of the system (29), when stabilising linear periodic state 



feedback of the form (27) is applied, is given by 

?fc+l = Mk + -BCT(fc)^'"0|fc = + B„(^f,)K^(^k))Ck 

or in general 

Cfc+i+l = ^a{k+i)S,k+i, i = 0,l,... 



where 



^a{k+i) = A + i?CT(fc+i)-^(T(fc+i) 



Then 



^k+jm+l 




£,k+jm+2 




Ck+jm+m 





^a{k+l)^a{k) 



^CT(fc+m-l) • • • ^a{k) 



Ck+jr 



and 



^^0|fc+jm+l 



^a(k)^k+jm 
K(T{k+l)S,k+jm+l 

-l^a{k+m—l)Ck+jm+in—l 
Ka(k) 

Ka(k+l)^a{k) 



Ka{k+m-l)^a{k+m-2) ■ ■ ■ ^<T{k) 



ik+jr 



(k) 



Thus, 



i=0 
oo m 



j=0 i=l 



|2 

m+i 1 1 q 



j=0 

oo 
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where 



Q = diag{q, g, • • • ,q) 

*<T(fe) = *a(fc+m-l) • • • ^<T(fc+l)^<T(fc) 

and 

oo 
j=0 

is a convergent series since the controller Kcr(k) is stabilizing. Thus P^^a(k) can be 
computed by solving the following Lyapunov equation 

Similarly, the sum of the control increments can be computed as 

oo oo m— 1 oo m—l 

Ju,k = ^ ll^'"0|A;+illP = ^ ^ ll^'"0|A;+jm+illP = ^ ^ \\^a{k+i)^k+jm+i\\f 
i=0 j=0 i=0 j=0 i=0 

oo 



— ^k+im^I(k]P'^crik)^k+jm 

oo 

j=0 

~ •^kPu,<T{k)-^k 

where 

R = diag(f, r, - ■ ■ , f) 



and 



Pu,a(k) = T.(Kik)f[Klik)RK,ik)]Ki 



r{k) 
i=0 

which can be computed by solving the Lyapunov equation 

Pu,a(k) = + K^(k)P'^a(k) 

Finally, let 

P^+u,a{k) = P^,<7{k) + Pu,a{k) 

oo oo 

= E(Kik)fi^lik)Q^^ik)]Kik) + 

j=0 i=0 

and it is clear that P^+u,a{k) can be computed as the solution of the Lyapunov 
equation 

Pi+u,a{k} = ^l{k)Pi+u,<Tik)^a{k) + ^l{k)Q^<T{k) + ^J(fe)^-^a(fc) 
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The cost (16), when the MMPC control law of the form (27) is applied, is given by 



(|30|) with 
q = 



' q 0' 


, f = 


" r " 











(35) 



and the corresponding initial conditions hold on ^fc; for example, at the beginning 
of an MMPC run it may be appropriate to set 



6 



Remark 4 Note that the cost defined in (30) differs in the first term from that 
defined in That is, in (30) we do not include any contribution from x^^, since 



that is fixed and cannot he influenced by the optimisation at time step k. 



Remark 5 It is seen that the system (29) is linear and periodic, while the cost (30) 
is quadratic with constant coefficients. Thus the optimal control law can be obtained 
from the theory given in and is of the form ([27]). This suggests yet another 
method to compute the terminal cost -^o-(fc)(2;A;+Ar|fc) to ensure nominal stability of 
MMPC in addition to that presented in Theorem [7J A family of MMPC designs 
may be obtained by optimising the cost function ( |30[ ) subject to the system (29) by 
choosing appropriate q and r matrices. 



Remark 6 // the second part of (^k> namely An^|^,_-|^, were included in the optimi- 
sation, so that the optimal cost became a function of Xk only, then the optimal cost, 
and the optimal solution, would be the same as that obtained with 'scheme 1 ' in our 
earlier papers, namely it would correspond to the cost resulting from allowing each 
'agent' to optimise all future inputs rather than just 'its own' input. 



5 Robust MMPC 

This section develops a robust version of MMPC. Uncertainty is introduced into the 
plant model as a bounded disturbance. The constraints which appear in the MMPC 
algorithm are then modified so that robust feasibility can be guaranteed, providing 
that it is achieved initially. 

The plant dynamics ([T]) are now extended to include an unmeasured but bounded 
disturbance Wk'- 

m 

Xk+i = Axk + ^ BjAuj^k + Ewk- (36) 
where Wk satisifies 

WkeW'ik (37) 
and W is a known, bounded set containing the origin. 
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As explained in Section [2| the system dynamics ( 36 ) can be re- written as a periodic 
linear system 

Xfc+i = Axk + -B^(fc)A% + Ewk (38) 

For robust MMPC we solve the following finite-time constrained linear periodic 
control problem, which we denote 'P„(^k){xk\k-, ^^\\k-v'^k~i)'- 



Minimize J*. = F^(^k){xk+N\k) + Eilo^ (Ikfc+iifcllg + \\uk+i\k\\l) 
wrt ^Uk+i\k {i = Q,m,2m,...,N-l) 

s.t. Auk+i\k ^l^i.aik), {i ^ 0,1, ■ . ■ ,N - 1) 

Xk+r\k ^ X,^a{k), {i = l,2,...,N ~1) (39) 

Xk+N\k € Tcj(k) 

Xk+i+i\k = Axk+i\k + -Bcr(fe+i) Aufc+j|fc, (i = 0, 1, . . .) 
^Uk+i\k = At2^.+,|fe_i + M^^„(k)Ewk^i, yi ^ jm 



Note some differences from Q. The predicted inputs and states are constrained to 
lie in sets Ui^cr{k) '^i,cr{k) which depend on how far into the prediction horizon they 
are, as well as on /c. The target set at the end of the horizon has been modified from 
Xx{K„(^k)) to Tfji^k)- Finally, the inputs which are not being optimised are assumed to 
be modified from their previously planned values by the feedback term Mi^„(^k)Ewk-i; 
note that the value of Ewk—i can be inferred from data {uj-i, Xj : j < k}. 

The constraint sets l^i^a{k)^ '^i,cr{k) ^-iid 7^(fc) are constructed to ensure robust feasi- 
bility, such that if some solution 



AU, 



k\k 



An*|^, An 



k+l\k' 



A-u 



fc+2|A:' 



Att 



k+N-l\k 



(40) 



is feasible at some time k then a candidate solution 

u{k+i)Ewk \ 



AUk, 



^^k+N-l\k~ 



■MN-2,a{k+l)EWk 



(41) 



\-Ka{k+l)X*kj^N\k+^'^N-l,cr(k+l)EwkJ 



is feasible at time k + \ for all -w^. G W. The designer chooses the feedback parame- 
ters Mi^„(^k) ^-iid offline (as in [27j, on which this development is based). 

To achieve this robust feasibility property, the state constraints Xfc G X are tightened 



usmg a recursion 



0,a{k) 



(42a) 



'^i+l,cr{k) — ^i,cr{k+l) ~ Li,a{k+l)EyV 



where 



-^0,cr(fc) — I 



(42b) 



(43a) 
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Li+i,a{k) = ^-^^i,<T(fc) + -B^(fc+j)^j,<T(fc) (43b) 

for the chosen feedback pohcy Mi^^i^f.-j and the "~" operator denotes the Pontryagin 
difference between two sets: 

A ^ B = {a \ a + b e A^b £ B} (44) 

Similarly, the input move constraint sets Au^ G are tightened using the following 
recursion 

^o,t7(fc) = ^aik) (45a) 



^i,aik) = Ui-l,a{k+l) ~ (45b) 

The terminal sets 7^(fc) have the robust invariance properties that, if x G 7^(fc) and 
w then 

(^ - Ba{k+N)Ka{^k+l)) X+ [AL^_i^^(fc+i) + -BCT(fc+Af)^Af-l,(T{A:+l)] Ew G 7^(fc+l) 

(46a) 



(46b) 



and 



(46c) 



The parameters Mi^^i^f.-j and ^C(j(fc) are chosen by the designer. The parameters ij,(T(fc), 
which relate the control perturbations in (41) to the corresponding changes in the 



state predictions, are then fixed by (43). These settings determine the amount of 



constraint tightening applied in (42). Typically, to achieve a large feasible region. 



the control policy chosen should minimise the quantities limited by the constraints. 

A restrictive but convenient choice of candidate policy is to select Mi^^f^f.-j, i = 
0, . . . , N — 2 such that ijAr,o-(fc) = VA; and then set Mjv_i,o-(fc) = 0, K^^^f.-^ = and 

%{k) = {0} yk. 

The following algorithm defines robust MMPC. It uses notations defined in ([T]), ([4]) 
and ( 23 ) . It is the same as Algorithm [T] except that problem ( 39 ) is solved instead 
of problem 



Algorithm 2 (Robust MMPC) 



1. Set k := ko- Initialise by solving problem (39), but optimising over all the 
variables Au|^_^_^f,,i = 0, 1, . . . , — 1. 

2. Apply control move Ano-(fc) = Au^j^ 
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3. Store planned moves A-u^ j^ifc. 

4. Pause for one time step, increment k, obtain new measurement Xk- 

5. Solve problem ( po] ). 

6. Go to step\^ 

We will need the following result concerning the use of the i^j^o-(fc) matrices. 
Lemma 1 Suppose that x^+i = + Ew^ and 

A{tfc+j|fc+i = Attfc+jifc + Mj_i^„(^k+i)Ewk, {j = 1,2,..., i) (47) 

Then 

Xk+i\k+i = Xk+i\k + Li^i^„{^kJ^i)Ewk, {i = l,2,...) (48) 

Proof: We prove the lemma by induction on i. Suppose the result is true for some 
i. Then 



Xk+i+l\k+l — ^Xk+i\k+l + B„(^k+i)^Ui.j^i\i.^i 



But 



Xk+i+l\k — + -Bcr(A;+j)^%+j|fc 

and, by assumption, 

^^f/c+ilfc+i = ^^k+i\k + ^i-i,^{k+i)Ewk 
so that 



(49) 
(50) 



(51) 



(52) 



Xk+i+i\k+i = + [ALj_i^^(fc_,_i) + i?^(fc+j)Mj_i^(^(fc+i)]i?w;fc (53) 

= Xk+i+i\k + Li^a{k+i)Ewk because of ( |43b[ ) (54) 



and hence the claimed result is true for i + 1. 

Now consider i = 1: Xf^^^j^^i = Xk+i = a^fc+iifc + Ewk, so the claimed result holds 



for i = I, since i>o,o-(fc+i) = by definition (43a 
Thus the result holds for z > 1. 



Theorem 3 If the system (36) is controlled using Algorithm^ and the initial opti- 
misation at time k = (ie step 1 of the algorithm) can be solved, and x kg G X, then 
(i) the optimisation remains feasible and (ii) the constraints £ X and AH^ G 



are satisfied for k > ko and for all disturbances satisfying (31). 
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Proof: (?) We will begin by showing that, by construction of the constraints in (42 ), 
feasibility at any time k implies feasibility at time k + 1. In particular, we will 
demonstrate feasibility by establishing that the candidate solution (41) satisfies all 
the constraints of the optimisation. Therefore, feasibility at time k = ko implies 
feasibility at all future times k > k^. 



Assume that we have a feasible solution (|40|) at time /c, and that An^|^ is applied 



as input to the plant (38). This results in the next plant state being 

Xk+i = Axk + B^(^k)^ul\k + ^^fe = Xk+i\k + Ewk 
Thus from (41) and Lemma [l] we have 



^fe+j+i|fc+i 

A^ffc+i+llfc+l 



Auk+i+i\k + Mi ,,^k+i)Ewk (i = 0, 1, . . . , iV - 1) 



(55) 



(56a) 
(56b) 



and since (40) was assumed feasible, we know Xk-^-i+i\k £ '^i+i,a{k) and Auk+i+i 



i+l,a{k) 



. Combining this with ( 56 ) , the definition of the Pontryagin difference ( 44 ) 



k e 



and the recursions (42b) and (45b), we know Xk+i+i\k+i € ^j,cr(fc+i) and An,t+i+i|fc+i ^ 



i,o-(fc+l) 



for all Wfc G W. 



We also need to show that Xk+N+i\k+i S 7^(fc+i) if the candidate solution (41) is 
applied. We have Xk+N\k ^ 7^(fc) by the assumption of feasibility at time k. 



Xk+N+l\k+l — Axk+N\k+l + B„(^k+N)^Uk+N\ 



But, from (41), 



k+l 



Ka(k+l)Xk+N\k + ^N~l,cj(k+l)EWk 



Attfc+Ar|A:+l = 

and, from Lemma [l] (since Au^+jvife = -K„(k+i)Xk+N\k 

Xk+N\k+1 = Xk+N\k + LN^i^cr{k+l)EWk 



(57) 



(58) 



(59) 



Hence, substituting (|58|) and (|59|) into (|57p gives 

Xk+N+l\k+l - 



[A - B^(^k+N)Ka{k+l)] ^k+N\k + 
+ [^-^Af-l,<T(fc+l) + B^(^k+N)MN-l,a{k+l)] EWk 



G Ta(^k+i) because of (46a 



(60) 
(61) 



Having established Xk+i+i\k e Xi+i,a(k)^ Auk+i+i\k e and Xk+N+i\k+i S 

Ta[k+i) for all Wk G W, the feasibility of the candidate solution has been proven, 
and thus the feasibility of the optimisation is proven. 

{ii) It remains to show that the state and input constraints are satisfied. Feasibility 
at all steps demands that Xk = Xk\k S X^ ^ji^ k) w hich from (42a) implies Xk G X. 
Similarly, Auk = Att^i^ E l^o,a{k) which from (45a) implies An^ G Ufc. 
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6 Examples 



This section demonstrates the potential benefits of MMPC by employing it in sim- 
ulation for the control of three different example systems. In the first example we 
consider nominal MMPC; we show how the cost formula can be used to evaluate 
some of the design choices. In the second and third examples comparisons are made 
between the robust MMPC scheme and standard — but also robustified — "syn- 
chronous" MPC (SMPC). In these examples all simulations were performed on the 
same PC with a 3.2GHz Intel Pentium 4 processor and 1GB RAM. Matlab ver- 
sion 7.1 (R14, Service Pack 3) was employed, using Simulink to simulate the system 
dynamics and the "quadprog" optimisation function to solve the necessary quadratic 
programming (QP) problems. Computation times were measured using the Matlab 
profiler. 



6.1 Nominal MMPC: Effects of Nu and updating sequence 

In this section, numerical examples will be given to illustrate how the cost formula 
for MMPC can be used for evaluating the effect of various values of Nu, and of the 
updating sequence, on the closed-loop performance, when constraints are not active. 



The cost formula for MMPC is calculated as (30) with (35) and initial condition 



of Cfc = [ 



xt 







Hence, only the upper- left n x n (n is the dimension of x^) 



sub-matrices are relevant in the cost computation. To be specific, the sub-matrices 
are P^^aik) and Pu,a{k) as shown below 



i+u,a{k) 



Pr 



a{k) 



P. 



x,a{k) 



u,a{k) 



u,a(k) 



■k 



where Pg+„ o-(/;), P^,cr{k) and Pu,a{k) are defined in (34), (32) and (33), respectively, 
and -k denotes a sub-matrix of compatible dimensions, which can be omitted from 
the cost computation. 

Therefore, the quadratic cost of MMPC can be computed as 



Ja{k) ~ Jx,a(k) ~^ 'Ju,a{k) ~ ^k Pu{k)^k 



where J^. 



^xMk) = xiPx,.^k)Xk and J„,.(,) = x^P„,.(,)X,. 

Now we have a way to compare different MMPC schemes, including differences in 
horizon lengths and update sequences. In other words, given two MMPC schemes, 
whose costs are Jj = x^PiXk and Jj = x]^PjXk, then 

Ji Jj — Xf. (fj Pj)xis 



Hence analysis of the properties of the difference Pi — Pj gives information on the 
relative merits of the two MMPC designs. For example. Pi — Pj > indicates that 
design j is better than design i for all initial conditions xq. 



20 



Consider the following two-input-two-output continuous-time plant 
1 1 



yi{s) 




_ y2{s) 





7s+l 3s+l 
2 1 

8s+l 4s+l 



ui{s) 

U2is) 



We chose the sampling time to be T = Is. For MMPC, the states were measured 
at T/m = 0.5s with ui and U2 alternatively applied at 0.5s intervals, but each held 
constant over a period of T = Is, that is, ui is updated at times (0, Is, 2s, . . .) and U2 
is updated at times (0.5s, 1.5s, . . .) (for the updating sequence of ni, U2, ui, ^2, • • • ,). 
For all the results listed below, the tuning parameters for MMPC are: q = I and 
r = 1. 

The horizon length is an important tuning parameter for MFC in general, and 
remains so for MMPC. With the MMPC cost formula, we can compute the cost and 
predict the performance difference with different N^. To illustrate, we generate a 
simulation scenario by adding a step input disturbance to the plant and see how the 
performance and cost vary with different A^„. This is done by replacing Am in ([T]) 
by {Au + Ad), namely modelling a change of input disturbance, so that an impulse 
on Ad corresponds to a step disturbance. This then allows such a disturbance to 
be represented by the initial condition x(0) = BAd, which in turn allows the use 
of formula ( |31| ) . (The results were checked against numerical estimation of the cost 
accumulated during simulation.) 

MMPC is a periodic control scheme; thus its performance depends on the time at 
which a disturbance occurs. More specifically, for the two-input plant considered 
here, depending on the time at which disturbances occur, ui^k may react first (ie 
control updating sequence {ui^k,U2.k,ui^k+i ■ ■ ■)), or U2^k may react first (ie control 
updating sequence {u2,k,ui,k,U2,k+i, ■ ■ ■))■ This depends on whether a{k) = 1 or 



cr(fe) = 2. This section uses the cost formula (31 ) to compare the cost of MMPC for 
these two different updating sequences in a specific scenario. 

Table [l] shows the eigenvalues of Pi — P2 as Nu varies, where Pi represents the cost 
matrix when a{k) = 1 (updating sequence (ni^^, n2,fc, . . .)) while P2 represents the 
cost matrix when a{k) = 2 (updating sequence {u2^k,ui^k, ■ ■ ■))■ From the table, 
it can be seen that Pi — P2 is indefinite, which means that one updating sequence 
is not definitely better than the other, but depends on the specific scenario — as 
expected. Fig. [2] compares the closed-loop performance between the two updating 
sequences when = 5. The solid lines show the response to a step disturbance on 
each input when ui is the first input to react to it (disturbance occurs at step k and 
cr(k) = 1), while the dashed lines show the response when U2 is the first input to 
react to the disturbance. The input trajectories approximately interchange in the 
two cases, as do the output trajectories, so there is little to choose between the two 
as regards performance. The cost difference xj^{Pi — P2)xk in this case is 0.3599, 
which means that the second updating sequence is slightly better than the first for 
this particular disturbance, as judged by the cost function. 



6.2 Robust MMPC: Spring-Mass Example 

This section considers the control of the simple mechanical system shown in Fig. [3| 
The system comprises four point masses moving in one dimension. Each has mass of 
five units and is connected to the adjacent masses by a spring of stiffness one unit. 



21 



Table 1: MMPC: Eigenvalues of (A - A) with different iV, 







Eigenvalues of (Pi — 


P2) 




1 


-9.4274 


-0.0069 


0.0000 


0.0000 


0.0042 


4.9050 


2 


-11.3972 


-0.0008 


-0.0002 


0.0001 


0.0014 


6.2401 


3 


-15.9939 


-0.0190 


-0.0001 


0.0001 


0.0289 


9.7714 


4 


-21.9017 


-0.0438 


-0.0002 


0.0001 


0.0603 


14.0665 


5 


-28.3386 


-0.0710 


-0.0001 


0.0000 


0.0905 


18.2002 



0.6 
0.4 
0.2 


-0.2 

















\ 






\ 
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time(s) 




10 20 30 

time(s) 




10 20 
time(s) 



Figure 2: Effects of updating sequence, step input disturbance, Nu = 5, solid: 
{ui,k,U2,k, ■ ■ ■), dashed: {u2,k,ui^k, ■ ■ •)• 
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y 




Figure 3: Spring-Mass Example System 



Each controller minimizes control energy subject to a constraint on the the position 
of mass 1, shown as output y in Fig. [sj Control energy is taken as / u{t)^ u{t)dt 
over a 400s simulation. The inputs are the control moves An^ applied to forces Ui 
acting on each mass, and therefore the control force levels u{t) are elements in an 
augmented state vector. All controllers were made robust to a disturbance force of 
up to 0.01 unit acting on mass 4. In the simulations, a disturbance pulse was applied 
to that mass of magnitude 0.01 from 50s to 200s. 

In the MMPC simulations, control moves were applied at intervals of one second, 
i.e. channel 1 moved at t = ti seconds, then channel 2 at t = ti + 1 seconds, and 
so on. In the comparison SMPC simulations, moves were made on all channels 
every four seconds, but to ensure fair comparison, the constraints were enforced at 
intervals of one second as in MMPC. Computation time is taken as the time spent 
in the "quadprog" function, totalled over all calls during the simulation. 

Figure |4] shows the control input signals and the output signals for each of the two 
controllers considered, using a horizon of 120s in both cases. The asynchronous 
control moves can be seen in the control signal plots from the MMPC simulation. 
In both cases, the output signal runs tightly against the constraint (shown dashed) 
for the duration of the disturbance pulse. This is as expected, since the objective 
is to minimize control energy and therefore the controller makes use of all available 
flexibility in the output constraint. The output under MMPC is slightly further 
from the limit than under SMPC, possibly because that controller effectively solves 
a more constrained problem due to the reduced decision variable set. However, the 
effect is not significant. 

To further illustrate the ability of the new robust MMPC to satisfy hard constraints 
despite disturbances, the simulation using MMPC was repeated using different con- 
straint levels. The resulting output signals are shown in Figure [5] In every case, 
the signal goes right to its limit, but never beyond, and the optimisations remain 
feasible. These results illustrate that the constraints are active in these simulations 
and that the robust MMPC method does not introduce undue conservatism. 

Table [2] compares detailed statistics from the results in Fig. [4| Observe that the 
performance, in terms of the control energy, is roughly the same for both controllers. 
However, MMPC is slightly faster than SMPC, since its sub-problems have only a 
quarter as many decision variables as SMPC. This illustrates the underlying premise 
of MMPC: it is faster to solve a sequence of four problems of 31 variables than one 
problem of 124. 

To further explore the issue of scalability, the simulations from Fig. |4] using SMPC 
and MMPC were repeated with various horizon lengths. Figure [6] shows the variation 
of total computation time with horizon length for both controllers. With a very short 
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□ 50 100 ISO 200 250 300 350 400 50 100 150 200 250 300 350 400 

Time Time 

(c) MMPC - Controls (d) MMPC - Output 

Figure 4: Spring-Mass Example: Responses to Disturbance Pulse 



horizon, SMPC is faster than MMPC. We hypothesize that this is due to overheads 
in the QP solver, such as set-up time, which dominate the solution time for small 
problems and therefore penalise the more frequent optimisation calls of MMPC. 
However, as the horizon length increases, the computation time becomes dominated 
by the actual solution process and MMPC scales more favorably than SMPC. 



6.3 Robust MMPC: Flight Dynamics Example 

This section considers longitudinal control of an A-7A Corsair II aircraft. The 
dynamics model was taken from Example 6.1 in Ref. [9] and augmented to in- 
clude a thrust input as well as the elevator input. Both inputs are constrained 
to [—0.04, 0.04] and the constraints are made robust to input disturbances in the 
range [—0.01, 0.01] on each channel. The simulation runs for 200s and a disturbance 
of 0.01 is applied to both channels from 20s to 120s. The planning horizon is 80s in 
all cases and the objective is to minimize where the state element X2 corresponds 
to the velocity normal to the aircraft axis in the body frame. 

Figure [7] shows the control and output signals from simulations using the two differ- 
ent controllers. SMPC executes moves on both channels at intervals of one second. 
MMPC performs a single move on alternating channels every half a second. Thus 
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Time 

Figure 5: Spring-Mass Example: Outputs from MMPC for Constraint Settings 0.2, 
0.4, 0.6, 0.8, and 1.0 

Table 2: Spring-Mass Example: Results for Each Controller Rejecting Disturbance 
Pulse. 



Controller 


SMPC 


MMPC 


/ u{t)'^u{t)dt X 1000 


4.312 


4.320 


Computation Time (s) 


6.6 


5.6 


N°- of QP Solutions 


100 


400 


N°' of Decision Vars. per QP 


124 


31 



the total number of moves on each channel in each simulation is the same. Table [3] 
compares the results using the same metrics as in the previous section, except for 
the performance which is here taken as the peak value of the normal velocity ||a;2||oo- 

Unlike in the spring-mass example, there is significant variation in performance 
between the two controllers. The MMPC controller, with its faster response time, 
is able to mitigate the short period response more effectively than SMPC, which 
leaves a significant spike at the onset of the disturbance, indicating that in this case, 
it is better to respond to a disturbance quickly with one channel than slowly with 
both. MMPC also requires significantly less computational effort than SMPC for this 
example. Note that the computation times are approximately in accordance with 
the expected 0(i/'^) behaviour, where u is the number of decision variables: in this 
example SMPC has 80 decision variables, and 200 QP problems are solved during 
the simulation, whereas MMPC has 41 variables, and 400 QP problems are solved. 
(200 X 80^) : (400 x 41^) = 3.7, which is quite close to the ratio of computation times 
42.25 : 9.15 = 4.6. 
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400 



Horizon (s) 

Figure 6: Spring-Mass Example: Variation of Computation Time with Horizon 
Length for SMPC and MMPC 



Table 3: Aircraft Example: Results for Each Controller Rejecting Disturbance Pulse. 



Controller 


SMPC 


MMPC 


\\X2(t)\\oo 


8.53 


6.49 


Computation Time 


42.25 


9.15 


N°- of QP Solutions 


200 


400 


N°- of Decision Vars. per QP 


80 
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7 Conclusion 

In this work a novel control scheme known as Multiplexed MPC was proposed, which 
is expected to be of practical benefit because it offers reduced computational com- 
plexity. Multiplexed model predictive control (MMPC) updates one input at a time, 
of a multi-input controlled plant. The motivation is to reduce the computational 
complexity of MPC, in order to allow reduced control update intervals. For some 
plants this leads to improved control, as a result of the controller being able to react 
to disturbances more quickly. MMPC scales well with increasing numbers of inputs, 
since the computational complexity depends only weakly on the number of inputs. 
The proposed MMPC scheme has been proved to be nominally stable. The nominal 
stability of a large class of other multiplexed MPC schemes follows by the same 
argument as we used in this paper. 

Some performance benefit over conventional MPC can be obtained as a result of 

faster reactions to disturbances, despite suboptimal solutions being obtained. This 
has been demonstrated by an example. However, the closed loop disturbance rejec- 
tion performance under MMPC is time varying because of the periodic nature of 
the control scheme. 
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(a) SMPC - Controls 



(b) SMPC - Output 



0.05 




20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200 

Time Time 



(c) MMPC - Controls (d) MMPC - Output 

Figure 7: Aircraft Example: Responses to Disturbance Pulse 



In this paper we have extended the basic MMPC idea to obtain robust feasibility and 
robust constraint satisfaction in the presence of unknown but bounded disturbances. 

Simulation examples have demonstrated that our scheme succeeds in maintaining 
constraint satisfaction and feasibility despite the presence of disturbances. Further- 
more, they have shown that performance improvements can indeed be obtained in 
some circumstances, compared with conventional MPC, they have indicated the kind 
of computational speed-up that can result from adoption of the MMPC scheme, and 
they have illustrated that these benefits are retained in circumstances where the 
constraints are active. 
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